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Abstract 

Differential elliptic flow spectra V2(pt) of n^, K^, p, A have been measured at yfs^N- 17.3 GeV 
around midrapidity by the CERN-CERES/NA45 experiment in mid-central Pb+Au collisions 
(10% of (Tgeo)- The pt range extends from about 0.1 GeV/c (0.55 GeV/c for A) to more than 
2 GeV/c. Protons below 0.4 GeV/c are directly identified by dE/dx. At higher pr, proton elliptic 
flow is derived as a constituent, besides n"^ and K^, of the elliptic flow of positive pion candidates. 
This retrieval requires additional inputs: (i) of the particle composition, and (ii) of V2{pt) of posi- 
tive pions. For (i), particle ratios obtained by NA49 are adapted to CERES conditions; for (ii), the 
measured V2{pt) of negative pions is substituted, assuming n"^ and elliptic flow magnitudes to 
be sufficiently close. The V2(pt) spectra are compared to ideal-hydrodynamics calculations. In 
synopsis of the series tt^ - K'^ - p - A, flow magnitudes are seen to fall with decreasing pj pro- 
gressively even below hydro calculations with early kinetic freeze-out {Tf- 160 MeV) leaving 
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not much time for hadronic evolution. The proton I'lipr) data show a downward swing towards 
low pt with excursions into negative V2 values. The pion-flow isospin asymmetry observed re- 
cently by STAR at RHIC, invaUdating in principle our working assumption, is found in its impact 
on proton flow bracketed from above by the direct proton flow data, and not to alter any of our 
conclusions. Results are discussed in perspective of recent viscous hydrodynamics studies which 
focus on late hadronic stages. 
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1. Introduction 



Among the prominent results from the Relativistic Heavy Ion Collider (RHIC) are observa- 
tions of strong elliptic flow in non-central collisions characterised by azimuthally 
anisotropic particle yields in the plane transverse to the beam direction (5, 6, Elliptic flow 
is quantified by V2, the second harmonic coefficient of the azimuthal particle distribution with 
respect to the reaction plane. The observations directly assert the importance of strong inter- 
actions among constituents of the expanding, hot and dense medium by which the geometrical 
anisotropy of the almond shaped overlap zone evolves into the momentum space anisotropy that 
is measured. This evolution is described by relativistic hydrodynamics (9). More specifically, 
the large V2 values agreed surprisingly well with predictions of hydrodynamics without dissipa- 
tion. This was interpreted as the early-time response of a locally equilibrated system of a very 
peculiar kind, the strongly interacting Quark Gluon Plasma (QGP), behaving as a nearly perfect 
liquid with a very small ratio rfj s of shear viscosity to entropy density (IQ, 11). 

From yfsNN- 200 GeV at RHIC to yfs^N- 2.76 TeV a bold step upward in nucleon-nucleon 
centre-of-mass energy was recently taken with the operation of the Large Hadron Collider (LHC). 
First measurements of elliptic flow in Pb-nPb collisions report viipj) for charged hadrons of sim- 
ilar magnitude and shape than at comparable centralities at RHIC ( T2, hlMj). The average V2 is 
si20% larger at the LHC but this increase is mainly due to the harder pj spectrum at LHC ener- 
gies. This is in agreement with hydrodynamic predictions extrapolated from RHIC data without 
change in the (very low) viscosity to entropy density dTsl as well as a hybrid calculation 
treating the QGP by ideal hydrodynamics and the late stages by a hadronic cascade model dlTb . 

At the Super Proton Synchrotron (SPS) energy, yfsNN- 17.3 GeV, elliptic flow magnitudes V2 
are about 30% lower flian at RHIC. The differential flow data V2(/?r) at SPS though 
strikingly similar in shape to the RHIC and LHC data, stay well below calculations of ideal 
hydrodynamics (^211 It should be noted that even at the highest RHIC energy some significant 
deviations remain (1221) . 

The failure of hydrodynamics at low energy, large impact parameters or large forward rapidi- 
ties has been ascribed to insufficient number densities at very early collision stages that hamper 
thermalization (l23l) . Strong dissipative effects are bound to set in after chemical freeze-out with 
growing mean free paths during the late hadronic expansion 

We present differential elliptic flow data V2(pr) of of strange particles A and for mid- 
central 158 AGeV Pb-i-Au colUsions collected by the CERES/NA45 experiment at the SPS Q; 
of negative pions complementing earlier, more peripheral data ( 19, 27); and of protons directly 
measured at low pj by dE/dx identification and retrieved at higher pr from the measured elliptic 
flow data of positive pion candidates containing besides pions and protons also kaons. The latter 
task assumes charge-independent pion flow, at least to an accuracy allowing to substitute the 
negative-pion V2(pt) measured for the positive-pion V2(pt) required. The particle composition 
also needed for the retrieval is fixed by particle ratios from the NA49 Collaboration that are 
adapted to CERES conditions. 

Elliptic flow data of identified particles as presented here are sparse at the SPS, especially for 
massive particles, and reaching down to low pr- We foresee these data to contribute valuable 
information on late stages of collective expansion which has been characterized by rescattering 
in the 'hadronic corona' (28) as an interplay between strong radial flow, elliptic flow, and thermal 
motion set by the freeze-out temperature (1211) . The proton flow data may turn out worth our effort 
as a probe specifically sensitive to hadronic viscosity (29). 

The assumed pion-flow isospin symmetry has recently been found invalidated by results of the 
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'beam energy scan' of the STAR Collaboration dsol) . in minimum-bias collisions at low ^JsMM 
and low pj- The yet unknown centrality dependence of the effect will be discussed in view 
of diverse physics scenarios. As a worst case remedy, our directly measured proton-flow data 
provide an upper limit for the uncertainty in proton V2 inflicted by the pion-flow asymmetry. 
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2. Experiment 



The 158 AGeV Pb+Au data were collected with the upgraded CERES/NA45 spectrometer 
during the heavy-ion run in 2000 at the CERN SPS. The CERES spectrometer is axially sym- 
metric around the beam direction and covers full azimuth at polar angles 7.7° < § < 14.7°, 
corresponding to a pseudorapidity range 2.05 < rj < 2.70 close to midrapidity (y,„id= 2.91); it is 
thus very well suited for elliptic flow studies. A cross-section view of the spectrometer is shown 
in Fig.[T| A detailed description of the CERES experiment is given in ('sT). 
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Figure 1 : The CERES/NA45 spectrometer during final data talcing in 2000. 



The radial-drift Time Projection Chamber (TPC) (32) is operated inside a magnetic field with 
maximum radial component of 0.5 T, providing a precise determination of the momentum. Par- 
ticle identification is achieved by the differential energy loss dE/dx along the tracks in the TPC. 
A doublet of radial Silicon Drift Detectors (SDD) (33), located at 10 and 13 cm downstream of 
a segmented Au target, is used for vertex reconstruction and tracking outside the field region. 
Charged particles emitted from the target are reconstructed by matching track segments in the 
SDD and in the TPC using a momentum-dependent matching window. The two Ring Imaging 
Cherenkov Counters (RlCHl, R1CH2) for electron identification were used in a previous CERES 
flow study of identified charged pions (Il9h . but not in the measurement reported on here. 



2.1. Momentum Resolution 

The momentum is measured by determining the deflection of the tracks within the TPC. The 
momentum resolution is therefore depending on the spatial track resolution, but is degraded by 
multiple scattering in the detector material. 

The results of an extensive Monte-Carlo study of the detector response were shown (l34h to be 
well approximated by the simple expression 



Ap/p = ^(2.0%)2 + (1.0% ■ /7[GeV/c])2. 



(1) 



2.2. Trigger Samples 

A sample of 30- 10^ events of 158 AGeV Pb+Au collisions was collected by a mixed-trigger se- 
lection with average centrality cr/cTgeo- 5.5%; this choice was made to enhance e+e" production, 
CERES' main objective. The track-multiplicity distribution for 'all triggers', shown in Fig.|2]by 
squares, has an average track number <Ntrack>u,- 157.9. It strongly deviates at low multiplicities 
from the minimum-bias distribution labeled {a) in Fig. |2] for the same average multiplicity, the 
minimum-bias distribution would have to be cut at Ntmck- 129.1 or at the top 11.4% of CTgeo- 

The limited statistics of our strange-particle spectra allowed for only two centrality classes. 
The 'top-central' part matches the minimum-bias distribution almost to its cut off at N,rack- 159 
and comprises with <Ntrack>\,„,,- 176.9 the top central 2.4% of (jgeo- The remainder are 'mid- 
central' triggers, <Ntrack>\,„ij- 136.2; these were used almost exclusively to collect the elliptic 
flow data. A precise definition of these triggers can only be provided by the distribution itself, 
together with the minimum-bias curve. For comparison to other experiments, we quote the sUce 
cut from the minimum-bias distribution that has identical average multiplicity: it extends from 
5.3% to 14.5% of CTgeo with a weighted average {cr/cTgeo} = 9.8%. 

2.3. Pion and Proton Identification 

Tracks in the TPC were reconstructed in the pseudorapidity range 2.05< 77 < 2.70. They had 
to pass quality cuts that required 

• transverse momentum pj above 50 MeV/c; 

• a minimum of 8 to 12 hits per track, depending on polar angle 




Figure 2: TPC track density for various trigger selections, (a) minimum-bias, (b) peripheral, (c) central. The combination 
of all triggers ('all triggers'), dominated by (c), has a weighted mean (cr/cr^^) of 5.5 %. It is split into 'mid-central' 
triggers (left, Ntrack < 159) of 9.8 % by which most of the flow data were collected; and 'top-central' triggers (right, 
Ntrack S 159 ) of 2.4 %, respectively. The horizontal scale ((r/a-ggo) on top applies to (a) only. See text. 
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Figure 3: Contour plots of the specific energy loss of charged particles vs momentum times charge sign. Full lines show 
Bethe-Bloch energy loss for n, K, p, dashed lines the ±1.5(T cut selecting pions. Relative dE/dx resolution, depending 
on number of hits per track, is typically 10%. 



• TPC and SDD tracks of charged pions and primary protons to match within a +3cr pj- 
dependent window. 

Pion candidates are identified by their specific-energy loss sampled along the tracks in the 
TPC. On average, there are more than 10 hits per track. 

In the two-dimensional scatter plot of Fig. |3] the measured specific energy loss dE/dx is 
shown as function of particle momentum p for both negative and positive charges. Shown are 
the cuts of Eq.|2] selecting pion candidates in a +1.5cr, i.e. ±15%, window (dashed lines) around 
the nominal energy loss for charged pions according to the Bethe-Bloch (BB) formula (full lines). 
The dE/dx cut is defined as 

0.85^(p,7r±) |bb< ^ip) lmeasured< l.l5^{p,n^) |bb ■ (2) 

dx dx dx 

Also shown are the Bethe-Bloch lines for kaons and protons, and it is obvious that over ex- 
tended ranges in momentum pions will be mixed with kaons and protons, in case of positive 
charge. Antiprotons are only about 6% of protons at mid-rapidity. At very low momenta protons 
are well identified by dE/dx as can be observed from Fig.|4] 

2.4. Reconstruction of A and K'^ 

The A particles are reconstructed via the decay channel A ^ p + with branching ratio 
BR = 63.9% and mean decay length ct= 7.89 cm (.35,) . Particle identification is performed 
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Figure 4: The low momentum part of Fig.|3]for positive charges with lines indicating the cut to select protons. 



using the dEjdx samples from the TPC by applying a ±1.5cr and +1.0cr window around the 
momentum dependent Bethe-Bloch values for pions and protons, respectively. The long decay 
length allows to accept only A decays which have no match to a SDD track within +3cr window. 
On the pair level, a pj dependent opening angle cut d^p,,* is applied. In addition, an Armenteros- 
Podolanski cut with qj <0.125 GeV/c and 0.0< a < 0.65 is applied in order to suppress 
background, admittedly with a considerable loss of signal. The a variable is a measure of the 
longitudinal momentum asymmetry, a - {ql—q])l{q^+qi), where q"^ and^^ denote longitudinal 
momentum components of and calculated with respect to p\ - + p^ . The qj variable 
is the momentum component of p^ in the transverse plane perpendicular to /?a. In the case of 
the A (/T^) particle one should exchange pf^ with pj^ (Pk°) in the above definitions. Fig.|5]is the 
2-dimensional a- qj scatter plot which shows the signatures of A, A and /T" reconstructed from 



the experimental data (l37h . 



2.4.7. A 

The combinatorial background is determined by rotating proton candidate tracks around the 
beam axis and constructing the invariant mass distribution. To decrease statistical errors in back- 
ground assessment, ten rotations by random angles are performed. 

The raw mass spectrum is shown in the left panel of Fig. |6] The pure A signal after subtraction 
of the combinatorial background shown at the right side has a non-Gaussian shape. This is owed 
to the facts that the observed mass and width depend on pj- andy and that the displaced secondary 
vertex is not used for recalculation of angles. The analysis is done separately in pj -y windows 
sufficiently small in size to keep the reconstructed A mass and width practically constant. The 
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Figure 5: Armenteros-Podolanski plot shows loci of A, A and reconstructed from experimental data; taken from ( BTI) . 
See text. 

signal distributions are fit by a Gaussian on constant backgrouncOl The mass of A particles 
strongly depends on pr but practically not on rapidity, while the width depends on both (l26|). 
Once the mass and width are established for a given y and pr both are kept constant for the rest 
of the analysis. 

With these cuts, values of the signal-to-background ratio S /B and the significance S / Vb of 
about 0.04 and 500, respectively, are obtained. Here, S stands for the signal and B for the 
background. Both quantities strongly depend on pr of the A. The largest values of S / yfB reside 
at pt ~ 1-1.5 GeV/c with y ~ 2; this is the most populated area in the y - pr plane of the 
reconstructed A. 

An example of reconstructed A in a given y-pT-(f> bin is shown in Fig. |7] (left). The yield of the 
A in a given bin is obtained by fitting the invariant mass distribution with a Gaussian. Plotting 
the yield versus for different pr and y values one can construct the dA^A/d0 distribution (Fig.|7] 
right). Fitting these distributions with a function c[l + 2v^, cos(20)], the observed elliptic flow 
values for different pr and y were extracted. The obtained v'^ coefficients were corrected for 
the event plane resolution as described in Sect. |3] 



2.4.2. 



The K'^ particles are reconstructed via the decay channel K'-^ — > n'^ + with branching 
BR = 68.95% and decay length ct - 2.68 cm (l35l) . In order to increase statistics, the dE/dx 
window is opened up to ±2 cr around the nominal Bethe-Bloch energy loss value for pions. 



found to be compatible with zero. 
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Figure 6: Left: A small enhancement of the signal is visible in the region of the A mass. Right: The invariant mass 
distribution of the A signal after subtraction of the normalized combinatorial background is significantly non-Gaussian. 
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Figure 7: Left: A reconstructed for 1.62 < y < 1.69 in rapidity, 0.675 < pr < 0.8 GeV/c in transverse momentum, and 
15° < (p < 30° in azimuth. Right: Elliptic flow pattern reconstructed from the A yield in bins for pj ~ 2.7 GeV/c. 



As the particle comes from the primary vertex, a possibihty to suppress fake track com- 
binations is given by a cut on the radial distance between the point where the back-extrapolated 
momentum vector of the K'^^ candidate intersects the x - y plane in the primary vertex. The nu- 
merical value of this cut is 0.02 cm. An opening-angle cut > 50 mrad is applied on pair 
candidates. Additionally, a cut on the z position of the secondary vertex (z > 1 cm) was applied. 
In order to suppress the contamination of A and A particles, an Armenteros-Podolanski cut with 
qr > 0.12 GeV/c was applied (see Fig.|5]). 

For subtraction of the combinatorial background the mixed-event technique is used. To pre- 
serve the event topology, only events with similar multiplicity and orientation of the event plane 
are allowed for mixing. Windows are set to ±10% and ±22°, respectively. The event mixing is 
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Figure 8: Left: The invariant mass distribution of tlie signal (red line) and the normalized combinatorial background 
(black line). Right: The invariant mass distribution of the K'i signal after subtraction of the normalized background. 



repeated 10 times. 

The K'-^ signal on normalized combinatorial background and after background subtraction is 
shown in Fig. [8] Mass and width of the reconstructed exhibit pr and rapidity dependences. 
Values of S /B and S / Vb of » 0.92 and 500 are obtained, respectively. 
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Figure 9: Left: K° reconstructed for 2.075 <y< 2.15, 0.2 < pr < 0.35 GeV/c and 45° < (f) < 60°. Right: Elliptic flow 
pattern reconstructed from the yield in (/> bins for pr ~ 1.7 GeV/c. 



In Fig. |9]an example of /T" reconstructed in a given y-pT-(p bin (left) and a K^^ flow pattern 
(right) is shown. The evaluation of the elliptic flow magnitude is done in the same way as for the 
A particles. 
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3. Elliptic Flow Analysis 

The flow analysis uses the event plane (EP) method, see e.g. (Q. BIS US) and further refer- 
ences therein. We give here only a short outline to clarify notations. The elliptic flow parameter 
V2 is the second term in the Fourier decomposition of azimuthal particle distributions in the plane 
transverse to the beam and with respect to the orientation of the reaction plane. However, the 
orientation of the reaction plane, being not known a priori, has to be reconstructed as 'event 
plane' (EP) for each event, 

. = A[ 1 + 2v^ cos(2(0,- - <^EP)) ]. (3) 
d((l>i - Obp) 

Here, 0, represents the azimuthal angles of outgoing particles. The anisotropy parameter vi, is 
smaller than V2 in magnitude due to the finite EP resolution. 

The azimuthal acceptance was divided into 100 adjacent slices a, b, c, d, a, b, c, d, a, ... 
such that every fourth slice was assigned to a subevent a, b, c, or d, respectively. To avoid 
autocorrelations, particle tracks employed for reconstruction of the EP and of were taken from 
non-adjacent slices only. 

Together with the reconstruction of the EP one calculates its resolution as the average differ- 
ence < Ofl - Oi > between the EP's reconstructed from two subevents a, b. Its inverse is the 
correction factor 'K given by 

7C = <2 cos(2(<l), - Oi))>-i/2 (4) 

by which the measured second harmonic is upcorrected, 

vi^-K v'^. (5) 

As the EP resolution depends on multiplicity, TC was calculated for different centralities. 
Fig. [To] shows for pion data the growing dispersion in EP orientation with multiplicity reflect- 
ing the fact that the decrease in anisotropy wins over the gain in statistics in deteriorating the 
resolution. 
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Figure 10: The correction factor vs TPC multiplicity for the 2-subevents method (pion data). 
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Figure 1 1 : Correction factors vs TPC multiplicity for A (full circles) and (open circles) flow analysis. 



To reduce autocorrelation effects, tracks chosen to be candidates for daughter particles were 
excluded from the determination of the event plane in case of A and . The event plane reso- 
lution was calculated for each multiplicity bin separately. The correction factors TC (Eq. |4|, the 
inverse of the resolution, are plotted in Fig. [TT]over TPC multiplicity for A'" and A. 

Due to instrumental inhomogeneities, the reconstructed event-plane density < dN/d<i>Ep > is 
not flat as it should be. In order to make it flat it is enough to apply successively first the method 
of recentering, and then the Fourier method of flattening 

(0). 



3.1. Correcting for HBT effects 

As we deal with a majority of charged pions, quantum effects among identical bosons give rise 
to space-momentum correlations of Hanbury-Brown and Twiss type (HBT): pions of the same 
charge tend to cluster in azimuth if their momentum difference - P2 1 is comparable to or below 
the uncertainty limit H/R. The critical momentum difference, with a typical source radius R x 
5 fm, is about 40 MeV/c. Since the average pion momentum is much larger, the effect is of short 
range in azimuth. The HBT effect correlates pairs of low relative momentum; it is a positive 
correlation which fakes genuine flow. 

In subtracting the HBT contributions to V2 we followed Ref. (1381) and used the source param- 
eters (in the standard Bertsch-Pratt parametrization) obtained from CERES HBT data 
The correlation function has been modified to take into account the effects of the Coulomb repul- 
sion (14 ih . Several iterations of the correction procedure were necessary to stabilize on the final 
value of the integrated elliptic flow; the latter decreased thereby in relative size by a; 10%. The 
results are shown in the next section. 

Concerning systematic uncertainties of the HBT correction, large relative corrections, and un- 
certainties are met with small magnitudes of flow at low pr', and while V2 quickly increases above 
0.5 GeV/c, the corrections diminish in relative size even more rapidly and so the uncertainties. 
An error estimate is reached by an educated guess based on a former study (Ref.28) where errors 
in the source parameters have been included. The chaoticity parameter was given a large error 
margin of ±50% in view of unknown influences of long-lived resonances, the momentum reso- 
lution and of pairs in which one or both pions are not primaries or of rho-decay origin. This lead 
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to an estimate of the relative uncertainty in the correction of +25%. The resulting relative errors 
in the HBT-corrected values amount to ±18% and +13% at pT= 0.25 GeV/c and 0.325 GeV/c, 
respectively. At pT= 0.50 GeV/c, the systematic error is down to ±3%. 
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4. Postanalysis of Elliptic Flow of Candidate Pions 



4.1. Overview 

As illustrated in Sect j2.3l our dEjdx cuts do not effectively filter out kaons and protons. In 
order to purify the elliptic flow data of tt" candidates (denoted "tt "), the knowledge of the K 
fraction as a function of pj is required, but also of K differential elliptic flow. We sketch the 
recovery of negative-pion elliptic flow in Sect J4.2l 

Positive pions are mixed with positive kaons and protons for p > 1.2 GeV/c. In order to isolate 
the proton elliptic flow from the measured V2 data of candidates, V2 , we will use the particle 
ratios at 158 AGeV which recently became available from measurements of Pb+Pb collisions by 
the NA49 Collaboration (42.) . 

The composition of the particle flux as accepted by the CERES spectrometer and filtered by 
the previous analysis cuts (26) has to be reconstructed. It is mandatory to properly simulate the 
effects of the pion-tuned dE/dx cut on kaons and protons. The dE/dx filter requires knowledge 
of particle momenta p- pr/sin 1?. Unfortunately, at the time this analysis was started, the infor- 
mation on polar angle § from the doublet of Silicon-Drift Detectors was no longer accessible. 
Without it, the resulting spread in sinj? over the acceptance, of almost a factor of two, would have 
blurred the dE/dx resolution. To avoid such degradation in quality, recourse was taken to a full 
Monte-Carlo (MC) simulation. We give details on the reconstruction of proton V2 in Sect j4.3l 

In a first approach to the acceptance correction we used analytic methods; although these 
calculations were made obsolete by the MC simulation, they serve as a valuable check on the 
final result. Besides, some are quite instructive and have been instrumental for preparing the 
input parameters of the MC simulations. 

4.2. Processing elliptic flow ofn^ candidates 

The expression linking the tt" differential elliptic flow Vj to the available data is 

vf = v'f + {v'f - vf ). (6) 

Here, is the flow parameter of K , r^ denotes the particle ratio r^ - NK-/Njr-. All quan- 
tities are functions of pj. Statistical errors enter in a way that forbids to apply standard error 
propagation. Therefore, a simulation was performed treating the three experimental inputs as 
Gaussian random variables with widths equal to their statistical errors. For every channel of the 
Pt spectrum to be incremented, the flow parameter V2{ptY is given by the average over many 
trials, its statistical error by the dispersion. 

As input for we use the differential elliptic flow data for which will be presented in 
Sect|7] Because of the very similar mass and quark contents of kaons, a possible difference in V2 
can be considered small compared to the present accuracy (see below). 

To improve on the statistical significance of the K^^ elliptic flow data, we have combined 
our results with those of NA49 '44). The latter data have been collected at the top 13% 
of CTgeo- The two data sets are shown in Fig. [12] together with a 1 -parameter fit according to 
viipr) - A Pj exp(-/:>7-). The best fit resulted in x^/ndf- 0.44. Shown is the best-fit curve for 
A = 7.18 ■ 10"^ in the centre, sandwiched between the ± Icr statistical error bands of ± 17% 
relative. 

To quantify the unknown admixture of K , we have used the K and tt" transverse momentum 



spectra measured by the NA49 collaboration (14211451) . Details are given in the following section 



for the analogous case of positive-pion candidates. All steps and procedures, like centrality 
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Figure 12: diflferential elliptic flow combining CERES data (this work; filled symbols) with NA49 data (ref.j^); 
open symbols). Centralities are (5.3 - 13)% of o-^,,, for CERES (average 9.8%) and top 13% for NA49 data. The three 
lines show the best-fit within the 1-cr errors bands. Errors are statistical. 



matching, dE/dx cut and acceptance corrections within the Monte-Carlo simulation, fully apply 
also to the negative-pion sample. 

4.3. Processing elliptic flow of candidates 

4.3.1. Outline 

To determine the differential elliptic flow of the minor component of protons from measured 
viipr) data of n'^ candidates, we make the simplifying assumption ~ ' ™d use the differ- 
ential flow of 71^ derived in Sect l4.2l as substitute for V2(n^) in the proton flow analysis below. 
For a discussion of possible violations of this assumption we refer to Sect l6.1.n and for our 
assessment of related uncertainties to Sect l6.2.l1 

The measured elliptic flow v"* of the ";7r^" candidate sample contains a contribution of proton 
elliptic flow 

v'f" = {N,. v\ + Nk* yf + Np vl)l{N„. + A^^^. + A^^). (7) 
More explicitly, the unknown magnitude Vj of the proton elliptic flow is derived as 

x'l = ((1 + rK. + rp) v'f - - r^. vf )lrp. (8) 

On the rh.s., we substitute the measured v'2(/r^) for v^iK^) ■ The only quantities yet unknown 
are the particle ratios = Nk* IN^* and - Np/N„t specifying the contents of and protons 
in the ";r^" sample, respectively. 

4.3.2. Diflferential particle yields 

Invariant yields of charged pions, charged kaons and protons at midrapidity for inelastic Pb-Pb 
collisions at yfsNN = 17.3 GeV, 

d^N(pT) ^ J_ d^N 
dydpT 2npT dydpr ,,^0 ' 
15 



10 



-6 



0.5 1 1.5 2 2.5 3 3.5 4 

p, (GeV/c) 



Figure 13: Invariant p7--diflferential yields of , , and protons from Ref. l l42h with statistical errors, shown 
together with 3-parameter fits obtained by minimizing data and curve are divided by 10. proton data by 100. The 

n" spectrum, very close to that of ;r^, is left out in order not to overload the figure. 



which are relevant for the present study are those displayed in Fig. 8 of Ref. The pt spectra 
are adjusted to CERES centrality as described below. For the purpose of sampling pj values at 
random which obey the proper density distributions, spectra are fit by 3-parameter exponential 
functions. The fits to the data are shown in Fig. [13] The authors of Ref. (42) have estimated 
the systematic errors, compounded from dEjdx, feed-down yields and acceptance corrections, 
as 2.2% for charged pions and K , 4.5% for , and 3.7% for protons. The fits also shown in 
Fig. [13] deviate form the data points by typically less than 1%. 

4.3.3. Matching Centralities 

Our mid-central collision window (5.3-14.5)% of cTgeo does not find a close match among 
NA49 centrality classes. Since data for the closest (5-12.5)% selection were not available ("4^, 
we have used a linear combination of centrality classes (0-5)% and (12.5-23.5)%. By visual 
inspection of the NA49 pt spectra for different centrality classes, a composition with equal 
weights seemed most appropriate. We also calculated the weighted means {cr/o-geo) for the two 
slices representing the NA49 centrality classes 1 and 3 to determine the required composition 
quantitatively: the mixture of 55% of class-1 centrality (0-5)% combined with 45% of class- 
3 centrality (12.5-23.5)% well reproduces the mean centrality of 9.8% of CERES mid-central 
triggers and was used. 
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Figure 14: The boundaries of CERES acceptance defined by polar angles i?= 7.7° and 14.7° shift away from central 
rapidity (>'„„rf = 2.91) with decreasing pj, the stronger the more massive the particle species, as shown here for pions, 
kaons, and protons. The //-acceptance 2.05 <ri < 2.70 coincides with y-acceptance at infinite momentum (fi= 1). 



4.3.4. Acceptance 

CERES accepts a cone in polar angle of 7.7° < i? < 14.7°. The small acceptance in pseudo- 
rapidity, Aj] = 0.65 units, spans the range 2.05 < tj < 2.70 and is close to midrapidity, ymid - 
2.91. For ptc smaller than particle mass mc^, the acceptance is shifted down in rapidity with 
decreasing pj as shown in Fig. [14] 

The practicable assumption is made that the doubly-differential yields of Eq. |9] factorize in 
y and pj- Particle yields within the CERES acceptance are then obtained by integrating the 
corresponding rapidity distributions between the /^^-dependent rapidity comers denoting the ac- 
ceptance for given m and pt, i.e. those shown in Fig. [14] We use the NA49 parametrization of the 
rapidity distributions by two identical Gaussians of width cr which are shifted from mid-rapidity 
by equal and opposite amounts +ys (47), 



ay 



(10) 



The values used in the present calculation for parameters N, cr, and are listed in Table [T] 
Those for charged pions and kaons were taken from Table III of Ref. (l47t) . The yield parameter 
for n'^ was downscaled relative to that for n'' by the same factor by which both the total multiplic- 
ities and the values of (dA^/dy)|y|<o.6 are observed to scale, by inspection of Table 11 of Ref. (1470 . 
For protons, the dNjAy parameters were determined by best-reproducing the distribution labeled 
CC2 in Fig. 6.6 of Ref. (1481) . This centrality class corresponds to (5-12.5)% of cr^^o, close to 
CERES centrality. 

Particle invariant yields have to be averaged over the acceptance before taking ratios; those 
will no longer reflect yields at mid-rapidity alone, but will also depend on the shapes of the 
dNI&y distributions with appreciable differences among the particle species. In addition, the shift 
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Figure 15: Acceptance correction functions A^cc of EQ. fTTl transforming NA49 particle ratios to CERES acceptance. 



in rapidity is taken into account. 

A suitable reference for the acceptance corrections of the invariant particle yields is the geo- 
metrical acceptance for massless particles (J3= 1) for which pseudo-rapidity coincides with rapid- 
ity, T] = y; Fig.[T4]would consist of two horizontal straight lines at y = -0.86 and aty - -0.210. 
With such reference, the acceptance correction is quantified by the ratio 



^acciPT, m) 



f 



yi'lhigi,; m,pT) jiy 

MAy 

y(ri,,„\ m,pT) "3' 



/ 



i-T-)dy 

= -0.86 ay 



(11) 



^we quote centre-of-mass rapidities downshifted from lab rapidities by the rapidity of the centre of mass in the lab, 
2.91. 





TT 


7r+ 


K- 


K+ 


P 




107.6±1.8 


104.4±1.8 


12.8±0.3 


23.4±0.6 


28.9±1.5 




1.18±0.02 


same as tt^ 


0.81±0.02 


0.88±0.04 


1.01±0.02 


ys 


0.72±0.02 


same as tt^ 


0.727±0.010 


0.839±0.012 


1.41±0.012 



Table 1: Parameters used to simulate dN/dy distributions using Eq. llQI From Ref. 1.47,) . See text. 
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The resulting correction functions are shown in Fig. [15] for the three particle ratios. The 
rapidity shift has opposite effects on the yield of protons and kaons: due to the rather compact 
distribution of kaons compared to the well-separated double humps in the dN/dy distribution 
of protons, the particle composition around 0.5 GeV/c is reduced by about 6% in kaons but 
enriched by 15% in protons over the /? = 1 reference. 

4.3.5. Monte-Carlo simulation 

Transverse momenta for given particle type are generated from the respective probability densi- 
ties, using standard methods i4% . Before doing so, pj spectra are adjusted to CERES centrality 
as described. 'Events' are filtered by conditions of acceptance and dE/dx cut. Once pr is cho- 
sen, boundaries of the acceptance for the respective particle species are defined in rapidity, as 
shown in Fig. [14] Then y is chosen at random from the density distribution of Eq. [TO]using the 
parameters of Table [1] Rapidity values outside the acceptance window result in rejection of the 
event. At this stage, the entire event is defined: the particle momentum p is calculated from 
m, pt, and y, which also fixes the Bethe-Bloch most-probable dE/dx value. Gaussian noise is 
added to simulate the experimental resolution cr= 0.10 dE/dx\BB- The historical filter set to (0.85- 
l.l5)dE/dx\BB of pions during data analysis is activated. Now, the survival fractions of kaons and 
protons are determined as error integrals over the dE/dx distributions between the /^^-dependent 
boundaries of the cuts on charged pions. The pr spectra are incremented by survival fractions /, 
0</< 1. 

Monte-Carlo generated spectra require subsequent normalization to the experimental data after 
which they were modelleqj. The NA49 particle yields are summed up from 0.30 GeV/c, the 
lower boundary of the first entry in the published spectra, to a reasonable upper boundary, i.e. to 
4.0 GeV/c. The normalization is achieved by setting the initial MC-yield, prior to any filtering 
and summed-up over the said pj range, equal to the sum over the respective data spectrum. 

The pr-integrated yield is multiplied by a factor xs to transform the mid-rapidity yield to the 
yield averaged over acceptance, 

xs = <dA^/d3;)f'^ES / aN/dy\y=o ■ (12) 

The transformation factors for particle ratios are listed in Table [2] Using identical binning, the 
statistical data errors are adopted. 



Since only particle ratios enter into Eq.[8] there is allowance for one free parameter common to all spectra. 



particle ratio 


K+/7r+ 


K-/7T- 


p/n^ 


Xs ratio 


1.401 ±0.026 


1.473 ±0.015 


1.415 ±0.016 



Table 2: Factors transforming MC-generated yield ratios at midrapidity to ratios of yields averaged over CERES' y- 
acceptance (Eg. 1 121 . Errors are propagated from Table[Q 
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5. Elliptic flow of charged pions: Results 



We present here the elliptic flow results of negative pions, illustrating the eff^ect of the HBT 
correction and the subtraction of the K component as outlined in Sect. 13.11 

The elliptic flow of negative pion candidates is shown in Fig. [T6l before and after the HBT 
correction. It is quite satisfactory to see that data start very close to zero and follow a quadratic 
Pt dependence for small pj- 
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Figure 16: Elliptic flow V2(pr) of rr candidates before (open circles) and after (filled circles) correction for the HBT 
effect (mid-central trigger). 

The /?7--differential elliptic flow V2{pt) of negative pion candidates is displayed in Fig. [l71 
together with the V2{pt) spectrum of negative pions obtained by subtracting the K component. 
Both spectra are corrected for HBT correlations as outlined in the previous section. Statistical 
errors are obtained by Monte-Carlo sampling along Eq. |6] treating all input data as Gaussian 
random variables. The estimate of 25% for the relative systematic error in the HBT correction 
amounts to about 15% of the corrected value around pT= 0.30 GeV/c (for details see Sect. 13. lb . 

The K impurity is seen to have only a minor effect. The correction starts out positive, raises 
the V2 value by about 0.005 at 0.6 GeV/c, decreases in size and changes sign at pr ~ 1.4 GeV/c; 
it turns negative where V2 of pions and kaons cross each other. The smallness of the correction 
follows from the rather small A^^- /N„- ratio and the two elliptic flow parameters being close in 
magnitude. 
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Figure 17: Diiferential elliptic flow V2(pt) of negative pion candidates (open circles); and of negative pions after remov- 
ing the admixture (filled circles). The vi spectra are corrected for HBT correlations. Statistical errors are compounded 
from all sources entering Eq. |6](mid-central triggers). 
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6. Elliptic flow of protons 



6.1. Systematic uncertainties 

Proton V2 parameters in the present study are obtained as differences between two 'large num- 
bers', the two pion flow measurements entering Eq. |8] Alarming as this might be in anticipation 
of relative errors getting out of control, the following evaluation will dispel such concern. A de- 
cisive advantage is that both measurements were performed and analyzed in identical settings^ 
In addition, CERES' full azimuthal coverage minimizes uncertainties in the EP determination 
and the flow measurement itself. 

Other uncertainties common to pions of both charges are related to the centrality determina- 
tion, and the dE/dx and acceptance cuts. Summing up the individual error estimates in quadra- 
ture yields 6%. A 'common-mode' error of this size is applied to both pion flow parameters and 
the response is calculated for the full pr range of interest. The deviations in proton V2 rise from 
negligible levels with increasing pr to Av2 ~ 0.008 at 2.5 GeV/c. 

Uncertainties linked to HBT correlations and the corrective measures taken are common to 
both pion flows and treated as relative errors. The HBT correction is largest at the lowest 
Pt bin of 0.325 GeV/c. A systematic relative uncertainty of 15% in the corrected values at 
Pt ~ 0.30 GeV/c causes an uncertainty in proton flow of Av2 < 0.004 which decreases quickly 
with rising pr to practically vanish around 0.6 GeV/c. The directly identified protons (the four 
asterisk points in Fig.[T9]below) are free of HBT correlations and related uncertainties. 

One may think also of uncertainties in only one of the two pion flow parameters which remain 
unbalanced. The subtraction of the K component from the ";r" candidate flow parameter is the 
only such example we are aware of. An estimated relative error of 8% in V2 of K combined with 
a 5% relative uncertainty in the K^/tt^ particle ratio has little impact. 

We have not corrected the pion elliptic flow data for non-flow correlations apart from HBT. 
Autocorrelation effects of short range between the samples used for flow measurement and for 
determining the event plane, respectively, were avoided by accepting only non-contiguous com- 
binations. In the proton flow analysis we may assume that on average non-flow correlations in the 
positive and the negative pion samples cancel each other; for the remaining proton elliptic flow 
we argue that jet-like correlations, the main physics source of non-flow correlations for pions, 
have negligible proton content at this energy. 

So much to systematic uncertainties in CERES data. The MC simulation introduces additional 
uncertainties. Errors in particle ratios have been derived from dn/dpr spectra of Ref.(42) as 5.0% 
for jn^, and 4.3% for p^ jn* . The MC-generated spectra carry small normalization errors, 
quoted in Table |2] Further errors arise by simulating dE/dx cuts and acceptance corrections, 
to our estimate of 3% each; they add up to 6.8% and 6.1% for the ln+ and pl7i+ ratios, 
respectively. The calculated uncertainties induced in proton V2 are largest at low pj but remain 
below 0.003. More harmful is the relative uncertainty in V2 of K^^: although deviations in proton 
V2 start well below 0.001, they approach 0.008 at the largest pT- 

6.1.1. A digression: pion flow asymmetry 

Our derivation of proton V2(pt) rested on the plausible assumption that V2{7r'^) and V2(7r ) 
are sufficiently close. With recent dicussions of exciting new physics related to chiral magnetic 



An exception is the different sign of traclc curvature in the TPC. The zero-deflection was fine-tuned by reconstructing 
the same mass for A and A )32l) . 



effects an asymmetry in elliptic flow of particles and antiparticles has come into focus 
and its observation was recently reported by the STAR Collaboration (,30i) . Preliminary data for 
minimum-bias Au+Au collisions at -iJsnn- 7.7 and 1 1.5 GeV show sizeable asymmetric^ 

^^(v2(n-)-V2(n^))/v2in^) (13) 

at very low 

There are reasons to assume the pion flow asymmetries wifl be considerably reduced towards 
more central collisions. We will go along some of the arguments that have been suggested, mostly 
by the authors of Refs. dsoi l30i) with emphasis on the expected pj and centrality dependence. 
We start with the more conventional scenarios. 

Resonance decays have been discussed for effects on pion elliptic flow (51) and possible vi- 
olations of the naive quark coalescence model (.52.) . Other than p-mesons, A(1232)-resonances 
may induce charge-dependent effects in flow by two reasons: (i) the isospin asymmetry between 
u and d flavour^ gives an edge to excitation of A° and A" over A^ and A^^, causing a surplus 
of tt" over tt^ decays, (ii) As A's are likely to be recombined from two hadrons (Nn), their 
V2 will increase over strict number-of-constituent-quark (NCQ) scaling, and decay pions are to 



show larger V2 than directly produced pions (l5lLl53h . Since A-decay pions are very soft, the flow 
asymmetry should show at very low pr- 

The asymmetry is transmitted only by decay pions which leave the fireball without being 
rescattered prior to thermal freeze-out. The condition is met by only a small fraction of A's 
decaying close to thermal freeze-outJ| The number of unscathed decay pions should be rather 
independent of centrality: production of A's grows with Nd, as do rescattering losses due to 
increased density. However, the asymmetry in the pion sample for more central collisions is 
reduced by an increasing share of thermal pions (oc Nch)- At the bottom line, ^ induced by 
resonance decays should scale inversely with < Nd, >■ For our 10% trigger we estimate a drop 
to about 40% relative to the (0-80)% minimum-bias trigger of the STAR data. 

Another scheme for violation of NCQ scaling has been proposed Jssl) in which partons carry 
larger amounts of flow strength arriving at midrapidity from stopped baryons, than those from 
qq pairs. This effect yields the correct sign, but is rather smallish: at yfsim- 8.86 GeV, the 
asymmetry is about 1%. At 17.3 GeV, it is further reduced by a factor 2.^. In this model, 
the fraction of constituent u quarks transported by baryon stopping decreases from 0.50 to 0.24 
between 8.86 GeV and 17.3 GeV, respectively!^ This suggests an equal reduction in proton 
density at midrapidity and hence of Coulomb repulsion of positive pions as a possible cause of 
the asymmetry, whatever its importance might be. 

The Chiral Magnetic Effect (CME): a transient magnetic field induces an electric current at 
finite baryo-chemical potential in presence of a chiral asymmetry between left- and right-handed 
quarks. The current generates an electric quadrupole with positively charged poles 'above' and 
'below' the reaction plane causing the asymmetry V2(;r") > V2(7[l) ( 1501) . The numerical estimates 
of (i5Q) are by and large consistent with the preliminary data (1301) . It is seen from Ref.(i56i) that 



^we changed the sign so that ^ > in accord with the STAR results. 

*we read from the J. Phys. paper /s^). Fig. 2 (top right) J?( values for Y*iviv= 7-7' 1 1.5, and 39 GeV and pT= 0.3 GeV 
of 30%, 26%, and 10%, respectively, which at pT= 0.5 GeV, reduce to to 11%, 8%, and 4%; minimum-bias tiigger 
^the neutron-to-proton ratio for Pb-Au is 1.52. 

^reduction of resonance yields relative to stable particles was described recently i Bih . 
'using the parameters of Ref. jSSi) and particle data from NA49 at 17.3 GeV. 
'"quoting values ot Xt defined in Ref.jS^. 
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Figure 18: Differential elliptic flow V2(pt) of n'^ candidates (open circles) and of identified negative pions (filled circles) 
shown already in Fig. [17] The two flow spectra difl'er by the and proton admixtures to the former, open-circled 
spectrum. A possible pion-flow asymmetry has been neglected. 



observables which measure the correlation of positive/negative charges to the reaction plane (15 7h 
are reduced by one order of magnitude between semi-peripheral collisions of 50-60% centrahty 
and mid-central collisions below 20% of the geometrical cross section^ It is plausible to assume 
that a CME-inflicted pion flow asymmetry is reduced also by one order of magnitude at our 10% 
centrality, compared to the minimum bias results. 

In summary, A-resonance decays and the CME seem plausible candidates for inducing an 
isospin asymmetry on pion flow at low pr, of the correct sign and possibly of the magnitude 
indicated by the preliminary STAR data (30) . The arguments presented suggest that ^ is con- 
siderably reduced in more central collisions compared to minimum-bias collisions studied by 
STAR. 



6.2. Results 

The primary data entering the proton flow analysis are displayed in Fig. [18] The indication 
of a negative excursion in V2 of the "tt^" candidates at low pr, absent in tt" flow, gave the 
incentive for the present reanalysis. The derivation of proton V2(pt) is performed by Monte- 
Carlo simulation of Eq.|8]using our viipT) data for K^^ and pj spectra of charged hadrons from the 
NA49 Collaboration. The simulation serves to adapt the NA49-based particle ratios to CERES 
conditions and not to introduce any model assumptions. 

The resulting proton V2{pt) spectrum between 0.30 GeV/c and 2.60 GeV/c is displayed in 
Fig. [19] It extends over the range where data both of particle spectra and kaon flow were avail- 
able. At lower pr, four additional points of directly identified proton^ are shown by asterisks. 



" refeiTing to the plotof a++,a i'.v centrality in Fig. 3ofRef./5^. 

'^see dE/dx plot Fig.|4]in Sect. [231 
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Figure 19: Proton V2(pt) (upper panel, filled circles) reconstructed from the "n^" candidate sample. The first four points 
apply to protons directly identified by dE/dx (solid asteiisks). Statistical errors derive from CERES V2("ff+") and V2(K°) 
data, and NA49 particle spectra, by Monte-Carlo sampHng. Mid-central events, (5.3-14.5)% of cr^gi,, weighted average 
9.8%. Lower panel: Systematic errors applicable to data points in upper panel, plotted to scale; see Sect. 16. II 



The first three of these data points, from pr - 0.05 GeV/c upward, are compatible with V2— 0, 
the fourth seems to bridge to the reconstructed points in an apparently abrupt downward swing. 
Considering the large statistical errors, the sawtooth impression may be misleading. The excur- 
sion at low transverse momenta of proton elliptic flow magnitudes below zero takes its minimum 
close to 0.4 GeV/c with V2= -0.0290+ 0.0092, 3.2cr below zeroH 

The systematic errors estimated in Sect J6.1l are displayed in the bottom panel of Fig. [19] The 
uncertainties due to the pion flow asymmetry are quantified in the following section lF.2.11 

6.2.7. Setting an upper bound on the flow asymmetry 

An asymmetry in pion elliptic flow threatens to falsify the results on proton elliptic flow we 
derived under the assumption that V2's are equal for n'^ and . Moreover, the subtraction of 
V2(7r") in place of the smaller V2{n^) in Eq. [8] would be an overcorrection that might have caused 
the peculiarity of the proton flow data, i.e. its turn to negative values at small pj. 

A meaningful upper bound on the asymmetry, however, is obtained from the proton V2 data: 
the lowest of the reconstructed points at pj- 0.325 GeV/c is close to the V2 point of directly 
identified protons at pj- 0.35 GeV/c, marked by an asterisk in Fig.[T9j being independent of the 
asymmetry issue, it serves as reference point. By tuning the asymmetry parameter in the Monte- 
Carlo simulation such as to move the reconstructed point by one standard deviation above the 



weighted mean of fourth direct and first reconstructed data point in Fig. 1191 
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Figure 20: Monte-Carlo simulation of proton fixing the asymmetry in pion flow J{ at 8.6% at pj= 0.325 GeV/c 
such that V2(0.325) (the first of the full-circle points) deviates by one cr from the directly identified point at 0.35 GeV/c 
(the last of the four asterisk points). At higher pj, the J{ values are made to decrease rapidly, modelling the strong pj 
dependence of the prehminary STAR data (full circles). The proton vi data of Fig. ll9l based on J{= are shown by open 
circles. See text. 



reference pointM. an upper boundJlmax{0.325)- 8.6% is obtained. To proceed further, we turn 
to the STAR data which show a dramatic decrease of ^ with pr, between 0.3 and 0.5 GeV/c by 
a factor of a; 3 (see footnote on p. 24). We use this pr dependence to extend our upper bound of 
asymmetry-related uncertainties beyond the calibration point: in the simulation shown in Fig.l20l 
y[ is taken to decrease linearly to 0.7% at pr- 0.525 GeV/c and stay constant above. We like 
these calculations be understood as defining the systematic errors of the proton viipr) data with 
respect to the pion-flow asymmetry issue. 

If the asymmetry were caused by A decays, the STAR data scaled to our centrality would come 
very close to our bound; and the strong pj dependence would be understood. If caused by CME, 
the scaled J[ would be much smaller than our upper bound. 



speaking about the rms of the combined data point enors 
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6.2.2. The rapidity window for proton flow 

Fig. |2T] demonstrates that protons for which the elHptic flow data \'i{pj) has been derived 
originate from about 1 .2 units below midrapidity. This reflects the rapidity shift of the acceptance 
for low-/?r protons while the r\ acceptance of CERES can be recognized in the tt^ distribution 
(reduced five fold). 
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7. Elliptic flow of A and K"^ particles 



7.1. Results 

The upper part of Fig. |22] shows the differential fiow data for Statistical errors are large, 
mostly due to the strong cut on the secondary vertex position during reconstruction. We show 
here the CERES data alone and refer to a comparison of the combined CERES and NA49 data 
with hydrodynamical calculations to Sect. |8] 

Absolute systematic errors in V2 of K'^, estimated by varying the cut on the z position of the 

secondary vertex, are !^ooo2 Pt ^ ^-^^ GeV/c and ^q q" for pr > 1.25 GeV/c. Systematic 
errors are quadratically small compared to statistical errors. 

The A elliptic flow displayed in the lower part of Fig.l22lshows a pj dependence characteristic 
for baryon elliptic flow. In the region of small pj, the magnitude of V2 is small but steadily 
increases with pj. At pr ~ 1.1 GeV/c, V2 exceeds 5% and rises further The absolute systematic 
error Av2 is estimated from two different ways of A reconstruction with emphasis either on the 
size of the signal 5, or on the signal-to-background ratio it is ^qqo' for pr < 1.6 GeV/c and 
for pt > 1.6 GeV/c which is again small compared to the statistical errors. 

For both species, /T^' and A, the lowest V2 value at pr = 0.175 GeV/c and 0.55 GeV/c, respec- 
tively, lies by about 1 cr below zero. 



7.2. Comparison to NA49 and STAR experiments 

A comparison to results from NA49 (44) at the same energy ( yfstm - 17.3 GeV) and to STAR 
results (15 8h at ^[sNN - 200 Ge is shown in Fig.|23] The NA49 and CERES data are in reasonably 
good agreement. In order to compare STAR to CERES results, the former have been rescaled 
to the centrality used in the CERES experiment. The appropriate factor is obtained by plotting 
the STAR V2 values vs centrality for different transverse momenta of A and K^^ particles. After 
rescaling, the STAR the V2 values measured at the RHIC energy are 15 - 20% higher due to the 
higher beam energy. 
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Figure 22: Differential elliptic flow V2(pt) of reconstructed mesons (top) and A hyperons (bottom). Mid-central 
collisions. Errors are purely statistical. 
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Figure 23: Comparison of A (top) and (bottom) elliptic flow data of CERES to measurements of NA49, Ref. f44h . 
and STAR, Ref. l58). 
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Figure 24: Pion differential elliptic flow compared to ideal hydrodynamics calculations ( 19, for two kinetic freeze- 
out temperatures, Tf= 120 MeV (blue solid) and Tf= 160 MeV (red dashed). Top: top-central collisions, Bottom: 
mid-central collisions. Data ai'e not corrected for admixture. Statistical errors. 



8. Comparison to ideal hydrodynamics calculations 

8.1. Overview 

We compare the results to ideal-hydrodynamics calculations by P. Huovinen These 



are performed in 2 -Hi dimensions assuming boost-invariant longitudinal flow. Initial conditions 
are fixed by reproducing the pr spectra of negatively charged particles and protons in Pb-i-Pb 
collisions at top SPS energy. Assumed is a first order phase transition to quark gluon plasma at 
a critical temperature of = 165 MeV. The calculations were done for two choices of kinetic 
freeze-out temperatureQ, T/ = 120 MeV and T/ = 160 MeV. 

Kinetic freeze-out at Tf- 160 MeV may be a handy way to shorten the evolution of the 
hadronic fireball and reduce vo thereby, but it fails in describing the proton spectra which 
come out too steep due to insuflicient radial flow, as noticed some time ago However, by 
comparing the data to alternative freeze-out conditions, we may find out how much the elliptic 
flow gains in magnitude during the hadronic fireball evolution, or even looses. 

8.2. elliptic flow 

The top-central data appear to be in perfect agreement for 'standard' T j- 120 MeV, as shown 
in the upper part of Fig. |24F^ However, it is hard to rule out, or it may even be likely, that the 



'also known as decoupling temperature Tdec 

'for easy comparison figures in Sect.|8]are plotted to the same scale in V2 and in pj. 
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Figure 25: K'^ elliptic flow data from Fig. [12] compared to ideal hydrodynamics calculations for Tf = 120 MeV (blue 
solid line) and Tf = 160 MeV (red dashed line). 



good agreement occurs by accidental cancellation between fluctuations that raise V2 and flow- 
damping effects reducing it. For mid-central collisions, the data seem to prefer a position in 
between the two hydro curves till about 1 .0 GeV/c and then saturate, whereas the hydro curves 
continue to rise about linearly with pr, see Fig.|24l bottom. Note, that until about 1.2 GeV/c, the 
mid-central pion flow data remain significantly above the Tf =160 MeV curve as to be expected. 

This is at variance with earlier CERES diff'erential flow data taken with more peripheral trig- 
gers {crla-geo- 13% - 26%) fl9). These elliptic flow data for and identified n"^, the latter with 
Pt threshold at 1.2 GeV/c, stay considerably below the predictions for T j - 120 MeV, but fall 
right on top of the Tf- 160 MeV lin^IZl In fact, the two CERES data sets seem to confirm that 
departures from ideal hydrodynamics increase with increasing impact parameter of the collision. 

Similar conflicts with ideal hydrodynamics have also been reported from RHIC ex pert - 
ments (f58[l60l) . Whether this indicates incomplete thermalization during prtmary stages (23l 6ll) . 
or increasing viscous corrections (|24[ 1291) . or a mixture of both, remains to be seen. 



8.3. K^^ elliptic flow 

The CERES/NA49 combined elliptic flow data for are compared in Fig |25] to ideal hy- 
drodynamics calculations. Surprisingly, for pj < 1.1 GeV/c data points tend to fall below the 
Tf = 160 MeV line. 

8.4. Proton elliptic flow 

The proton elliptic flow V2{pt) data are shown again in Fig. |26]for compartson with hydro- 
dynamics calculations. With standard freeze-out temperature, the calculation overpredicts the 
data as in the pion case. What appeared to be a tendency for K^^ is for proton V2{pt) plain fact: 
the majority of data points resides below the Tf = 160 MeV line. With decreasing pr, the ideal- 
hydro curve bends smoothly towards zero while the data continue to fall about linearly until 
about 0.2 GeV/c. For such early freeze-out, as implied by the large Tf, it is hard to attribute 



referring to Fig. Ic in JTgh . 
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Figure 26: Proton v^ipr ) data of Fig. [T9] compared to ideal hydrodynamics predictions for two freeze-out temperatures 
(see Fig.l24l. 

the V2 reduction to a decreased contribution from hadronic interactions. The observation rather 
seems to suggest that a suppression mechanism or dissipation is at work which is not present in 
the ideal-hydro calculation. 

To what extent negative V2 values are reached is surely a matter of the asymmetry value; for 
our upper-bound J[, all negative V2's with one exception remain negative. But more important, 
the perception that V2 data for protons towards low pr fall further below and away from the 
160 MeV-ideal hydro curve as in the case of K^^ remains fully valid. The theory seems to miss 
some basic ingredient, at least in the present modelling, to account for the marked reduction in 
V2 at low pt akeady indicated for K^^. 

8.5. A elliptic flow 

The reduction in V2 is seen in Fig. |22]to persist for A elliptic flow. The A results could be 
termed 'perfectly in line' with previous findings, would that not be an overstatement in view of 
the large statistical errors and the sparse low- pr data in case of the A. Even an excursion into 
negative V2 values as seen for protons would fit into the A flow data at very low pr- 

The data suggest that deviations from ideal hydrodynamics grow with particle mass. That in 
fact the deviations seen in A flow are not stronger than those in proton flow may well be due to 
the small lever arm: the relative gain from m(p) to m(A) is only 19%, compared to the steps m{7T) 
to m{K) (250%), and m{K) to m{p) (90%); and to limited data precision. 



Hydro, T =120 MeV 

Hydro, T =160 MeV 
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Figure 27: A elliptic flow data and predictions of ideal hydrodynamics for Tf= 120 MeV (blue solid) and Tf = 160 MeV 
(red dashed). Mid-central events. Errors purely statistical. 
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9. Discussion and Conclusion 

We have presented differential elliptic flow measurements vaCpr) at y/s^N^ 17.3 GeV of 
and A in mid-central collisions, supplemented by V2(f r) of negative pions. We presented also 
differential elliptic flow V2(pt) of protons, at low pj by direct identification, for pj >0.35 GeV/c 
by reconstruction from impure positive-pion samples. The synopsis of differential elliptic flow 
spectra from n to to p and A in comparison with calculations of ideal hydrodynamics discloses 
a marked decrease of V2 values towards low pj, getting stronger the larger the hadron masses. 
That negative values are reached in proton V2 maybe seen as the most prominent feature of this 
trend. 

These results were shown to be robust with respect to isospin asymmetries in elliptic flow of 
charged pions as reported recently by STAR at RHIC dsol) for minimum-bias Au + Au collisions. 
An upper limit on J[ could be derived by using our direct proton point for reference, at the lowest 
Pj measured where J[ is supposedly largest. 

Indications of negative-valued elliptic flow of heavy particles from SPS and RHIC experiments 



at low pt and central to mid-central collisions have been reported (621 I63LI64I) . but differential 
flow data at SPS or RHIC energies — with more than one or two points below zero — have to our 
knowledge not been presented before. A tendency of massive particles towards negative V2 values 
at low Pj is ascribed to the conjunction of elliptic flow and strong radial expansion under specific 
freeze-out conditions. The authors of Ref. (12 lb like to illustrate and quantify these phenomena 
using a 'blast-wave' description (65) which is based on an interplay of transverse flow, imparting 
a momentum gain proportional to particle mass, and elliptic flow anisotropy which renders such 
momentum shift larger in-plane then out-of plane. The particle density at low pj is depleted 
in-plane more than out-of plane, and elliptic flow values tend to become negative. The depletion 
mechanism works best for steep pj spectra, i.e. those with low inverse slopes or temperature. 



In a recent hydrodynamic study (l29h a large hadronic shear viscosity-to-entropy ratio was im- 
plemented to calculate pj spectra and differential elliptic flow of hadrons produced in 200 AGe V 
Au-hAu collisions. While viscosity in the late hadronic phase suppresses elliptic flow in general, 
the consequences are striking for protons: for rfj s > 0.42 and Tj- < 120 MeV, proton V2 turns 
negative at small pj and the similarity to our data can hardly be overseer[3. The authors stress 
that this effect is entirely owed to viscous corrections 6f (24) creating large average transverse 
pressure. The viscous mechanism, growing with nij, appears to initiate a kind of blast wave with 
the potential to produce effects very similar to those we observe in proton elliptic flow. 

Some caution is advised. An earlier but related viscous hydrodynamics study (66) remained in- 
conclusive as proton spectra and charged hadron elliptic flow were found to have non-compatible 
requirements on the size of (constant) rj/s. Besides, Ref.(67) points to possible consequences of 
yet unresolved ambiguities in viscous hydrodynamic calculations for identified particle observ- 
ables. 

We hope that our results are useful to better understand the role of shear viscosity during the 
late hadronic stages that terminate heavy-ion collisions at all energies. 



we refer to Fig. 3d and Fig. 4d. of Ref l l29b . 
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